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ABSTRACT 

We present the results of a long-term observation campaign of the extra- 
galactic wind-accreting black-hole X-ray binary LMC X-l, using the Propor- 
tional Counter Array on the Rossi X-Ray Timing Explorer (RXTE). The obser- 
vations show that LMC X-l's accretion disk exhibits an anomalous temperature- 
luminosity relation. We use deep archival RXTE observations to show that large 
movements across the temperature-luminosity space occupied by the system can 
take place on time scales as short as half an hour. These changes cannot be 
adequately explained by perturbations that propagate from the outer disk on a 
viscous timescale. We propose instead that the apparent disk variations reflect 
rapid fluctuations within the Compton up-scattering coronal material, which oc- 
cults the inner parts of the disk. The expected relationship between the observed 
disk luminosity and apparent disk temperature derived from the variable occulta- 
tion model is quantitatively shown to be in good agreement with the observations. 
Two other observations support this picture: an inverse correlation between the 
flux in the power-law spectral component and the fitted inner disk temperature, 
and a near-constant total photon flux, suggesting that the inner disk is not ejected 
when a lower temperature is observed. 

Subject headings: Accretion, X-rays: Binaries (LMC X-l, LMC X-3, Cygnus X-l), 
Black Hole Physics, Stars: Winds 
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Introduction 



LMC X-l is (with Cyg X-l) one of only two known persistently luminous x-ray binaries 
consisting of a black hole accreting the wind of a massive blue star. The 10.91 ± 1.41M 



black hole is in a 3.9 0917 ±0.00005 day orbit flOrosz et al 



( 1 Cowley et al. 



20091 ) about an 07 III companion 



19951 ) . As the companion both drives a strong wind and is far from filling 



its Roche lobe, wind accretion feeds the black hole 



which accretes at an average of O.I G Lp.m (IGou et al. 



low/hard state ( Wilms et al 



200lh 



Gou et al 



Nowak et al. 



20011 ). The system, 



20091 ) has never been observed in the 



( 120091 ) used x-ray data, including some of 



the RXTE observations we use here, to derive a spin parameter a = 0.92(+0.05, —0.07) for 
the black hole in LMC X-l using fits to the disk blackbody component of the spectrum. 

LMC X-l's persistent occupation of the high/soft state contrasts with the behavior 



of Cyg X-l , which it closely re sembles in other res pects. Cyg X 



black hole (IHerrerq et al. 



companion ((Bolton 



19951 ) in a 5.6 day orbit (jPaczynskil Il974f ) about an 09.7 lab 



harbors a ~ lOA-f, 



© 



19721 ). Though Cyg X -l is also a wind-accr eting; system, its typical 



luminosity in the high state is ~ 0.03LEdd (jGierlihski et al. 



1999), which is much lower than 



LMC X-l's. This implies that the black hole in LMC X-l accretes at a higher rate than 
its counterpart in Cyg X-l, perhaps because the companions' wind speeds and densities 
differ, or that accretion mechanisms with different radiative efficiencies are taking place. In 



addition, Cyg X -l exhibits we 



low. 



1975 



r ard state (lAgrawal et al. 



Sanford et a 



( iNowak et al. 



2001 



-documented transitions between the high/soft state and th e 



1972 



Baity et al. 



1973 ; 



Heise et al. 



1975 



Holt et al. 



1976 . 



1975). LMC X-l , in contrast, displays persistent soft-state emission 



Wilms et al. 



20011 ). While both LMC X-l and Cyg X-l's softest states 
can be fitted with a power law plus disk blackbody model, they also differ markedly. The 
traditional soft state model, in whic h the disk blackbody dominates the e nergy spectrum, 



accurately fits LMC X-l's spectra ( lEbisawa et al. 



1989 ; 



Yao et al. 



2005|). Cyg X-l's 
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spectra, however, are energetically dominated by the power law in both the hard and soft 
states. Thus, while these two systems have similar companions and orbital properties, they 
differ in their spectral characteristics and evolution. 

An even more interesting comparison can be made between LMC X-l and LMC X-3, 
a black-hole binary that accretes via Roche-lobe overflow at comparable luminosity, as a 
means of highlighting the observational differences between disk and wind accretion. Much 
of this paper will be devoted to pointing out and interpreting these differences. 

LMC X-3 is a persistently bright x-ray binary in the Large Magellanic Cloud, and is 
the only other system with a dynamically confirmed black hole. As such, it provides a useful 
counterpoint to the wind-accreting black hole binary systems mentioned above . LMC X-3 



is_a 9M black hole ( Ivan der Klis et al. 



19851 ) with an evolved B5 companion (ISoria et al. 



200 lh . A detailed listing of the system's dynamical properties, along with those of LMC X-l 
and Cyg X-l, is provided in Table [TJ Differences between the systems are expected to 
reflect differences in the accretion flows resulting from wind versus Roche-lobe overflow 
accretion: LMC X-3 accretes high angular momentum material through its first Lagrange 
point, while LMC X-l accretes low angular momentum material from its companion's 
stellar wind. LMC X-3 and LMC X-l share spectral signatures of a persistent accretion 



disk, thoug 



i LM C X-l's disk spectrum defies the L ex T 4 relation ( JMitsuda et al 



Dunn et al 



20101 ) 



1984; 



This paper examines how LMC X-l's spectral properties and evolution fit with, 
and possibly extend, existing black hole binary accretion models. In particular, our 
results ultimately indicate that LMC X-l's spectral behavior is consistent with sporadic 
obscuration of the innermost part of a stable accretion disk. 

Section 2 of this paper presents the data processing and spectral fitting techniques 
employed in obtaining the results presented in Section 3. Section 4 discusses how LMC X-l's 
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Table 1:: System parameters for LMC X-l, LMC X-3, and Cygnus X-l, with references 
in parenthesis. The Eddington luminosity for each source was calculated from black hole 
mass. 



Distance 

(kpc) 



Black Hole 

Mass (M ) 



System 
Inclination 

(degrees) 



Companion 
Type 



Companion 
Mass (M Q ) 



Companion 
Radius 

(Re) 



Orbital Pe- 
riod (days) 



Average 
Soft State 
Luminos- 
ity (% 
^km) 



01 



(2001) 



LMC X-l 



48.10 ±2.22 (1) 



10.91 ±1.41 (I) 



36.38 ±1.92° (1) 



07 III (8) 



31.79 ±3.48 (1) 



17.0 ±0.8 (I) 



3.90917 ±0.00005 (I) 



16% 



Orosz et all fl2009fk 



5(5) 

D; (9) 

!3|); (12) 



Herrero et al. 



Nowak et al. 



LMC X-3 



52.0 ±0.6 (2) 



11.1 ±1.4 (4) 



50.0° < i < 70.0° (12) 



B5 IV (6) 



4.0 < M < 4.7 (6) 



4.4 (6) 



1.70479 ± 0.00004 (5) 



< 30% (9) 



2) Idi Benedetto! (11997 ) 



19951): (6) 



Soria et al 



van der Klis et al 



toiLiio) 



(I1985I ) 



121 



Cyg X-l 



2.5 (3) 



10.1 (5) 



36° < i < 67° (11) 



09.7 lab (5) 



17.8 (5) 



22.7 ±2.3 (7) 



5.6 (3) 



3% (10) 



Paczvnskil (11974 ) : (4) 



teooih: m 



Gierlihski et al. 



Ziolkowski 



(120051): f8) 



(|1999|); (11) 



Gierl hski et al 



Gou et al. 



Davis fc Hartmann 
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unique spectral properties can be reconciled with existing accretion flow models, while 
Section 5 presents our conclusions and possible avenues for further investigation. 

2. Observations and Data Analysis 

The data analyzed in this work fall into two categories. The bulk of the observations 
come from our twice- weekly RXTE monitoring campaign of LMC X-l and LMC X-3. The 
LMC X-l data were complemented by long, individual RXTE observations available on 
the HEASARC. These data were reduced and analyzed using identical processes. The two 
classes of LMC X-l observations prove useful for probing different timescales of spectral 
variations. 

2.1. Observations 

The main observing campaign consisted of brief, uninterrupted observations of 
LMC X-l and LMC X-3, conducted twice a week since August 2007 and March 2006, 
respectively. These observations ranged in length between 506 and 5470 seconds, with 
average lengths of 1661 seconds for LMC X-l and 1877 seconds for LMC X-3. The 
LMC X-l observations were offset from the source by 15.764' to reduce contamination from 
the nearby pulsar PSR B0540-69. These twice-weekly observations give reliable pictures 
of the typical X-ray spectra of LMC X-l and LMC X-3 over significant periods of time. 
While this is a major strength of these data sets, it also a weakness, in that they only have 
the capacity to reveal spectral variations that emerge over a minimum of several days. 

To examine shorter timescale variations in LMC X-l, we take advantage of the 
seventeen deep archival RXTE observations listed in Table [2j We "chopped" these data sets 
into shorter spectra (less than 90 min) corresponding to one orbit of the RXTE spacecraft. 
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These chopped archival observations can reveal spectral changes that emerge over fractions 
of a day. 



2.2. Data Analysis 

All data were reduced using standard HEASOFT v6.7 tools. Data taken within 
30 minutes after the satellite exited the South Atlantic Anomaly (SAA), or when the 
source was at less than 10° elevation above Earth's limb, were excluded from subsequent 
analysis. All observations only considered data from the second proportional counter 
unit (PCU2) of RXTE's Proportional Counter Array (PCA), due to its consistent 
performance throughout the mission. The mission-long faint background model 
pca_bkgd_cmf aintl7_eMv20051128.mdl was used in this analysis, as the total count 
rate for both LMC sources lies below 40 photons/s/PCU. Data were extracted using the 
Standard2 data mode. 



We used uniform XSPE C v!2.5.1 ([Arnaud 



Per RXTE recommendation ( jJahoda et al. 



1996J ) fitting procedures for all data. 



20061 ). systematic errors of 1% were added 



in quadrature to the Poisson counting noise errors. Because the detector response is 
poorly understood within the first three PCA channels, we exclude those channels from 
our analysis, resulting in a lower energy threshold just under 3 keV. Above 16 keV, the 
signal to noise ratio for both sources is too low to provide useful fits to the data, so our 
fits also exclude all data for energies over 16 keV. A 29% decrease in the PCA efficiency 



(JJahoda et al 



20061 ) relative to an on-axis source resulted from the 15.764' off-axis pointing 
of our LMC X-l observations. All fluxes derived from the off-axis LMC X-lobservations 
were multiplied by a factor of 1.4 to account for this efficiency factor. Representative 



spectra, model fits, and fit residuals for LMC X-l and LMC X-3 are presented in Figures [T] 
and 121 respectively. 

Spectra were fitted with a disk blackbody plus power law emission model. The first 
spectral component is a multicolored blackbody spectrum generate d by the accretion d isk, 



which is accounted for usin g the standard dis kbb XSPEC model (IMitsuda et al. 



The simpl XSPEC model (ISteiner et al. 



1984h . 



20091 ) accounts for the power-law emission, which 



extends to high energies, and which is thought to be caus ed by Compton up- scattering 



19761). A pure 



of disk photons through a population of hotter electrons (IShapiro et al. 
power-law model runs the risk of promoting spectral fits that inaccurately trade disk 
luminosity for coronal energy flux. The XSPEC simpl model avoids this by turning over 
at low energies, and by convolving an arbitrary input photon spectrum with the specified 
Compton scattering prescription. It is, therefore, less prone to overestimating the amount 
of non-thermal emission. Simpl can account for either pure photon up-scattering, or for 
both up- and down-scattering. Because the average energy of coronal electrons is two orders 
of magnitude higher than the highest disk blackbody temperature, coronal emission from 
our sources can be safely approximated by pure up-scattering. 

Finally, to account for photoelectric absorption within the intervening ISM, we convolve 
the compound model w ith p habs. Our a p proac h does not conflict with the findings of 



Levine fc Corbet 



(120061 ) and 



Hanke et al. 



(120 101 ). which detect orbital phase variations in 
the absorption column towards LMC X-l. These orbital phase variations suggest that LMC 
X-l's companion drives a significant wind. While the variable column density significantly 
effects LMC X-l's observed soft x-ray emission, it does not impact the higher energy range 
examined in this analysis. 

Table [3] lists the fitting parameters and results from the following fitting procedure. 
The simpl power law index, the absorption column, and the inverse Compton scattering 
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Fig. 1. — : normalFit and data-to- model ratio for LMC X-l data fit by the multi-colored 
disk and simpl (power-law) spectral model. 
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Fig. 2. — : Fit and data-to-model ratio for LMC X-3 data fitted by multi-colored disk and 
simpl (power-law) spectral model. 
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flag are all frozen a t the listed values. The absorption column densities are drawn from 



NowaketaL 



( 1200 ll ) for both LMC X-l and LMC X-3. The innermost disk temperature 
starts, but is not frozen, at the physically reasonable energy of 1 keV. The spectral fits 
output the inverse Compton photon scattering fraction, the innermost disk temperature, 
the overall normalization of the diskbb model component, and \ 2 goodness-of-fit statistics. 
Once the model has been fitted to the data, it is used to calculate both comprehensive and 
component-specific photon and energy fluxes. 



The low photon counts from both LMC X-l and LMC X-3 limit the usable energy 
range in our 1.5 ksec observations to no higher than 16 keV. With counting statistics of this 
quality, it is not possible to simultaneously constrain both the power-law index and the 
inner-disk temperature. Since our observations cannot determine all six of the parameters in 
the phabs x simpl x diskbb model (see table), we freeze the power-law index, which is the 
most stable at long time scales. The deep archival observations of LMC X-l, prior to being 
chopped into orbit-by-orbit sections, have sufficient photons to simultaneously constrain all 
disk and power-law fit parameters. We use the average power law index derived from these 
long observations, which is consistent with a single value, to fix the power law index for 
all of the single-orbit LMC X-l spectra. For LMC X-3, we employ the average power law 



index found in the same way by 



Smith et al 



20071 ). In both cases, the values for all deep 



pointings are statistically consistent with the average. 



Results 



We find that LMC X-l remained in the soft state over the entire observation 
program, and has an anomalous disk-temperature- versus- luminosity relation relative to the 
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expectation of a disk blackbody. 

Figure [3] demonstrates how the luminosity of a typical accretion disk, such as the one 
present i n LMC X-3's soft st ate, follows the expected modified L ex Tf n Stefan-Boltzmann 



relation ( JMitsuda et al. 



19841) • The line in Figure [3j which represents this relation, was 
not fitted to the data beyond finding a suitable normalization constant. The relation 
nevertheless fits the data very well, aside from the highly uncertain points at low luminosities 
and high temperatures. These points come from observations taken after Feb. 23, 2009 
(MJD 54885), when LMC X-3 transitioned from a disk-dominated high/soft state to the 
low/hard state. The fitting routine we employ is not appropriate for the low/hard state, 
during which disk emission is suppressed and the power-law index hardens significantly. The 
poorly constrained values that appear in the lower right corner of Figure HI well away from 
the Stefan-Boltzmann relation, result from using a high/soft state model to fit low/hard 
state data. 

Figure H] similarly shows the disk temperature- luminosity relation for LMC X-l. The 
dashed lines are two normalizations of the same Stefan-Boltzmann relation used in Figure 
HJ and clearly cannot fit the data. 

Figure [5] compares how the inner disk temperatures of LMC X-3 and LMC X-l vary 
over time. The figure shows only a subset of the total observing campaign in order to better 
display short term variations from each source. The change in disk temperature from one 
observation to the next is markedly more continuous for LMC X-3 than for LMC X-l. 
Over longer time periods, however, LMC X-3's inner disk temperature evolves significantly. 
Such long-term disk temperature evolution is absent in LMC X-l. 

Figure O shows how LMC X-l's inner disk temperature and luminosity varied over the 
course of two of the archival observations listed in Table |2j The bottom panel shows data 
from the longest observation, while the top panel shows data from an observation with a 
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Fig. 3. — : The disk temperature-luminosity relationship observed in LMC X-3. The dashed 
line is the modified accretion disk Stefan-Boltzmann relationship. The points at low lumi- 
nosities and high temperatures, which lie far from the model, are due to a transition to the 
low/hard state towards the end of the observing campaign. 
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Fig. 4. — : The disk temperature-luminosity relationship observed in LMC X-1. The dashed 
lines are the modified accretion disk Stefan-Bolt zmann relationship shown in Figure [31 and 
differ only in their overall normalization constants. 
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typical number of exposure intervals. 

The inner disk appears to evolve over a period of several hours. Rather than evolving 
along the Stefan-Boltzmann relation, the disk temperature and luminosity fall along a line 
of a very different slope. The order of the points is not uniform along this line either, 
suggesting motion faster than has been resolved by the data. 

Figure [7| shows LMC X-l's variability in disk, coronal, and total luminosities, as 
well as in the total photon emission rate, over the course of our observing program. The 
mi ld discrepancy be tween LMC X-l's overall accretion luminosity as presented here and 



m 



GouetaL 



( 120091 ) comes from the narrower energy range we employ in calculating the 
source's overall flux. Both analyses, however, find similar disk-to-coronal flux ratios in LMC 
X-l. Percent variation, defined as standard deviation /mean, quantifies the extent to which 
the values in each panel vary. The system's disk luminosity shows 14.5 % variation, while 
its coronal up-scattering luminosity shows 51.8 % variation. The system's overall x-ray 
luminosity has 15.2 % variation. The total photon flux from LMC X-l is the system's most 
stable feature, with only 12.5 % variation. While the total number of photons emanating 
from the system remains relatively steady, the amount of energy that the corona contributes 
to the total emission varies significantly. 

Figure [8] displays the relation between LMC X-l's inner disk temperature and the 
inverse Compton scattering fraction. While observations with cooler disk temperatures and 
higher scattering fractions suffer larger fitting uncertainties, a Pearson correlation analysis 
confirms the existance of a statistically significant negative correlation between these two 
parameters, even when the two upper-leftmost points in Figure [8] are removed. Since the 
scattering fraction maps directly to the inverse Compton optical depth, Figure [8] shows that 
the inner edge of the disk appears cooler when there is more coronal material. 
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4. Discussion and Conclusions 

LMC X-l's apparent defiance of the Stefan-Boltzmann relation may stem from 
variations in the disk itself, or from variations in the scattering medium surrounding it. 
Because the total number of photons is more nearly constant than the disk luminosity (see 
Figure [7]), variations in the number of photons diverted to the power law by the corona 
seem more likely to explain the anomalous temperature/luminosity relations than variations 
within the disk itself. 

The first step in quantitatively distinguishing between these two explanations for the 
anomalous disk relation is to compare the observed disk variation timescales to timescales 
inherent to the disk and to the corona. In Figure 0, LMC X-l 's disk randomly samples 
a wide swath of the temperature parameter space in the space of a few hours. Individual 
points in the archival observation are separated by an average of 2320 s, which provides an 
estimate of the timescale for apparent disk variations. We can compare this value to the 
disk's viscous timescale to assess whether the disk's structure is capable of changing over 
such a brief interval of time. 



The disk's viscous timescale scales with both radius and viscosity, as (IFrank et al 

2OO2I) 



tvisc ~ (1) 

V 

The circularization radius, R C irc, of the accreted wind material gives a lower limit for the 
outer disk radius. The disk radius, and therefore the viscous timescale, may actually be 
larger due to outward angular momentum transfer. Since we are primarily concerned with 
finding a lower limit on the viscous timescale, however, R C i rc provides a conservative value 
of the disk radius. The circularization radius is defined as the radius at which the accreted 
matter's angular momentum due to Keplerian motion e quals the angular momentum it 



carried upon initial capture by the black hole. Following 



Frank et al. 



( 120021 ) we assume that 
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the wind's initial angular momentum about the black hole goes as 



(2) 



where R acc represents the black hole's gravitational capture radius, and u is the black hole's 
orbital angular velocity about the companion. From there, we obtain 



lL r j rr 



G 3 Ml 



BH 



UT 



(3) 



J rel 



where Mbh represents mass of the black hole, and 



2 _ 2 | 2 
V rel — V w + V orb- 



(4) 



Therefore, the disk's size is strongl y dependent on th e 



unde rgoes line-driven acceleration ( 



Castor et al. 



1975 



speed o 



Owocki 



the stellar wind, which 



1994: 



20001). The wind's ve locity increases with distance from the star (ILamers fc Casinelli 



Kudritzki fc Puis 



1999 



Kudritzki et al 



19891 ) as 



v w (r) = v c 



0.9983 



-R* 



(5) 



Here, i^ is the wind's terminal velocity and R* is the companion's radius. The factor of 
0.9983 is chosen such that v w equals the escape velocity at the star's surface. The values 



of Voo and /3 are poorly constrained f or O-star winds, but 



t ypically range from 



000-2000 



19991 ). We 



km/s and from 0.5-1.5, respectively (IKudritzki et al.lll989l ). (jLamers fc Casinelli! 
conservatively assume Voo = 1700 km/s, in order to compa re our results with other works 



that consider clumping in O-star winds (IDucci et al 



20091 ). We carry the effects of the full 



range of /3 values through all subsequent calculations to obtain a range of results. 

The final step in calculating R C i rc is determining v w where the wind reaches the 
accretion cylinder. To find R acc , we equate the wind's gravitational and kinetic energy 
terms with respect to the black hole and find 

2GM BH 



Rn 



(6) 



J rel 
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Solving (T4|) , (J3J) and ([6]) numerically over the full range of /3 values constrains v w to lie 
between 154 and 635 km/s by the time it reaches the black hole's accretion cylinder. R C i rc 
is thu s limited to values between 6.7 x 10 8 and 2.4 x 10 10 cm, in agreement with previous 



work ( Nowak et al. 



20011 ). 



Having found R aC c, we turn to calculating the disk's viscosity, v. The Shakura-Sunyaev 
a-disk prescription posits 



v = aCeH 



where c s and H are the disk's local sound speed and 



(7) 



height, respectiv ely. Observations 



20071 ). Like LMC X-l, 



of transient x-ray binary outbursts indicate a ~ 0.1 (jKing et al. 

x-ray transients while in outburst host ionized accretion disks around a central black hole. 

These similarities motivate us t o use a = 0.1 in the following calculations. Substituting 



equilibrium thin-disk formulae (IFrank et al. 
(EJ) gives 



20021 ) for the sound speed and disk height into 

3/10 



v = 1.8 x lO 1 ^ 5 ^ 1 V 1/4 # 



3/4 
10 



Rr, 



R 



1/2 



cm /g. 



(8) 



Here, mi is Mbh in units of M Q , M\e is the accretion rate in units of 10 16 grams per second, 
Rio is radius in units of 10 10 cm, and R m i n is radius of the innermost stable orbit. Simple 
geometry requires that the accretion rate depends on the companion's mass loss rate as 

2 



*» = i (£= 



^sys 



M*g/s 



(9) 



with a typical O-star M* of 10 -5 M Q per year. The factor of four in the denominator 
reflects the the fact that the black hole presents a circular cross-section to the spherically 
expanding wind. Equation ($$) shows that the accretion rate onto the black hole ranges 
between 2.2 x 10 17 and 9.6 x 10 18 g/s, depending on the chosen value of (3. These values are 
in good agreement with accretion rates derived from the observed x-ray luminosity using 



L ~ 0.1c 2 M BH , 



(10) 
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which range betwee n 8.4 x 10 17 and 2 .0 x 10 18 g/s. The 0.1 efficiency factor in (fTUj) 



is recommended by 



Frank et al 



(120021 ). Equation [TU1 arises from the condition that all 



gravitational potential energy up to the last stable orbit is radiated away. 

Using ([[]) to combine the results of ([S} and ([3]) gives 5 x 10 4 < t V i SC < 4 x 10 6 seconds. 
These values of t visc are at least an order of magnitude longer than the ~ 10 3 second 
timescale of inner disk variation indicated by Figure [6l This provides the first piece of 
evidence that the observed disk luminosity variations are not caused by a process inherent 
to the disk itself. But the calculation of the disk viscous timescale above assumes that 
the entire disk is involved in any change. The Roche-lobe accreting black-hole binary 
GRS 1915+105 shows very fast spectral c hanges in many of its variability states near the 



Eddington luminosity (jBelloni et al. 



20001 ). which appear to be repeated ejections of the 
innermost parts of the disk. But when the inner disk is ejected and the spectrum turns hard 
in GRS 1915+105, the total x-ray count rate also drops dramatically. It is the absence of 
any such change in count rate in LMC X-l that leads us to favor obscuration by marginally 
optically thick, ionized coronal material to an ejection mechanism. 

The slope of the temperature-luminosity relation that the chopped archival observations 
in Figure [6] occupy provides the second indication that LMC X-l's apparent disk luminosity 
variations are caused by a process unrelated to the accretion disk. If the variations in disk 
luminosity were due to changes in the accretion rate, the points would move sequentially 
along the Stefan-Boltzmann relation shown in Figures [3] and HJ 

Since the anomalous disk temperature-luminosity relation is poorly explained by 
changes within the accretion disk, we ask whether changes in the corona could explain the 
relation. First, we consider which coronal properties are compatible with the observations. 
Then, we speculate how such a corona might form. 

Figure [8] provides a clue as to how the coronal material factors in LMC X-l's anomalous 
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disk temperature-luminosity relation. The fraction of disk photons that suffer inverse 
Compton scattering while traversing the corona is lowest when the apparent inner disk 
temperature is around 0.9 keV. This value lies at the higher end of the observed range of 
disk temperatures for this source. The inverse Compton scattering fraction increases with 
decreasing inner disk temperature, up to a scattering fraction of nearly unity at Tj n of 
~ 0.4 keV. Such behavior is consistent with a steady disk whose inner radii are sometimes 
obscured by a cloud of energetic electrons. In that case, the number of disk photons, which 
form the seed photons for inverse Compton scattering in the overlying corona, would remain 
constant. The fraction of original photons that suffer up-scattering would depend only on 
the amount of coronal material present. Figure [7] shows that LMC X-l's total photon flux 
remains largely constant, while its coronal luminosity varies substantially. 

To quantify th is model, we integrate the emission from the a— disk flux prescription 



(IFrank et al. 



20021 ) outwards from inner radii at temperatures ranging between 0.4 and 
0.9 keV. Figure [9] shows the resulting temperature-luminosity relation, as fitted to all of 
the chopped archival datasets with four or more sub-intervals. The disk temperature and 
luminosity are constant in each observation. The amount of inner disk obscuration is the 
only variable that determines where the chopped observations fall along the best-fit line. 
The model appears consistent with all the data, particularly with observation 30087-01-04, 
which provides the strongest constraint. Note that these individual chopped observations 
run roughly perpendicular to the Stefan-Boltzmann relation in temperature-luminosity 
space. Suzaku or NuSTAR observations may allow the detection of variations of the 
power-law index on relevant (hour) timescales; hardness variations might then be interpreted 
as variations in Comptonizing optical depth, in which case hardness might be expected to 
correlate with the flux in the tail. 



Figures [9] and [8] raise questions about the analysis in 



Gou et al. 



(120091 ). which derives 
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the spin of LMC X-l's black hole from measurements of the inner disk temperature. 
Their analysis requires that the inner disk always be unobscured. Our analysis, however, 
shows that the degree of inner disk obscuration in LMC X-l varies significantly on short 
timescales. In systems such a s LMC X-3, whe re the disk is not obscured, there is no 



difficulty with the approach of 



GouetaL 



(bOQflh 



Why, in fact, does LMC X-3 have less apparent inner disk obscuration than LMC X-l 
does? We can imagine at least three possible answers. The first option is that the coronal 
optical depth is much larger, in general, in LMC X-l than in LMC X-3. This possibility is 
ruled out by the overlap between the two sources in terms of the flux ratio between the disk 
and power-law components; when the disk luminosity in LMC X-3is in decline, a significant 
power-law tail can form in that system without di srupting; the excellent Stefan-Boltzmann 



relation for the disk component (ISmith et al 



20071 ) . A second possible answer supposes that 



the geometries of the coronae are very different: LMC X-3 has a geometrically large corona, 
while LMC X-l has one that is more centrally concentrated, so that it obscures more of the 
inner disk while producing the same amount of upscattering. A third possibility supposes 
that the corona in both cases has a conelike shape (possibly a jet, but not necessarily, from 
our evidence alone), and that the apparent differences between the two sourc es are due to 



an in clination effect. In LMC X-l, which has an inclination of 36.38 ± 1.92° (jOrosz et al. 



20091 ). the conical corona may block the inner parts of the dis k as we look down throu gh it 



while in LMC X-3, with an inclination between 50° and 70° (Ivan der Klis et al 



19851 ). the 



cone presents itself in semi-profile, where it doesn't obscure the inner disk but we can still 
see the Compton upscattered flux it produces coming out its sides. 

While either of the latter two pictures is consistent with our data, the last is more 
appealing for two reasons. Firstly, it does not require an unexplained difference in the 
form of the corona in the two binaries. Secondly, it may someday be possible to verify 
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or to disprove using x-ray polarization observations. A 'down the barrel' jet observation 
(e.g. LMC X-l) would sport a lower over-all scattering polarization fraction than an 
observation of the same type of geometry viewed 'in profile' (e.g. LMC X-3) would. The 
precise magnitude of the expected polarization fraction difference will govern how easily 
its presence or absence can be verified. In any case, future observations with the GEMS 
x-ray polarimeter will at least begin to explore the relevant paramet er space, proy ided the 



instrument can spend upwards of 10 5 seconds observing each source ( IJahoda 



2010|) 



Finally, we turn to the question of how a population of energetic electrons forms and 
dissipates above the innermost region of the disk on timescales shorter than t visc . Figure 
[10] illustrates one possible accretion scenario that could account for LMC X-l's unusual 
properties. 

In this picture, the incoming wind material sh eds its net angular mom entum in the 



'post shock" before accreting onto the black hole ( jHoyle fc Lyttleton 



19391). Most of 



the shocked gas remains in the disk plane, and enters the accretion disk (middle arrow). 
A smaller portion of the heated post-shock matter launches on randomly distributed 
radial trajectories towards the black hole (upper and lower arrows). This component of 
the accretion flow rapidly heats on its radial descent into the black hole's gravitational 
potential well. Unlike the disk, which is dense enough to cool efficiently via bremsstrahlung, 
the diffuse infalling matter cannot efficiently cool via free-free interactions. By the time 
the radially accreting matter passes over the inner disk, it can have reached electron 



temperatures a 



Dove t 



( iNarayan fc Yi 



1995 



re 100 keV required for the obs e rved inverse Com pton upscattering 



Ichimaru 



1977 



Rees et al. 



1982 



EsinetaL 



19981 ) . This model allows 



the corona to change independently of the disk. 

The presence of two accretion flow components, which are largely independent beyond 
the post-shock, is consistent with the constant observed total photon number and minimal 
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disk luminosity variation in LMC X-l. Changes in the wind mass flux, due perhaps to 
wind clumping, register immediately in the radial, coronal flow. Viscous diffusion in the 
disk, meanwhile, smooths out those same sharp mass discontinuities. While the coronal 
surface density, and therefore the total energy of coronal emission changes rapidly, the disk 
emission evolves more slowly and less dramatically. The scenario of independent flows with 
the thin disk remaining intact beneath the corona, even when the latter is optically thick, 
is consistent with the picture derived from obser vations of the hard state and hard-to-soft 



20021 ). Recent observations of 
ativistic line profi les 



(Miller et al. 



2006). 



transitions in IE 1740.7-2942 and GRS 1758-258 flSmith et al 

the iron fluorescence line in black holes in the hard state have shown re 

indicating that the inner disk remains present even in the hard state 

until extremely low lumniosities, below 1% of Eddington, are reached ( iTomsick et al. 

In the luminosity range of LMC X-l, we therefore expect the inner disk to remain present 

whether the thin disk is actually continuous or wh ether its inner p ortion is recondensed 



2009|). 



from a disk that is disrupted at intermediate radii (JLiu et al 



20 lit ) 



Differing patterns of x-ray variability (e.g. Figure [5]) show promise as a means of 
identifying black hole companions when optical observations cannot. Fast changes in the 
inner disk temperature (Figure [5]), and an anomalous temperature/luminosity relation 
(Figure HJ), both combined with a stable net photon flux (Figure |7J), may signal the presence 
of a wind accretor, perhaps combined with a low inclination, just as long-term hysteresis in 
state changes appears to be a signature of Roche-lobe overflow accretion and a large disk 



( ISmith et al. 



2002 



20071). 



The type of two-component accretion flow proposed in this paper is best suited to black 
hole binary systems where the companion drives a high-velocity stellar wind. At the same 
time, this particular model leaves ample room for other coronal generation mechanisms that 
have already been proposed for black hole binary systems accreting via Roche-lobe overflow. 
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At least two types of observations can either help to confirm or to disprove the 
hypothesis we have proposed. The first focuses on recent x-ray observations of IC 10 X-l 



and NGC 300 X-l (Barnard et al 



20081 ). Both systems have disk to power-law luminosity 
ratios, as well as high-mass main-sequence companions, that resemble LMC X-l's. Binary 
systems that contain both an O-type star capable of driving such a wind and a black hole 
companion are rare, which has made this particular combination of accretion mechanisms 
particularly hard to observe. Further observations of IC 10 X-l and NGC 300 X-l, 
however, can potentially determine how deep their apparent similarities to LMC X-l 
run, and will provide more opportunities to test our model's veracity. The second test 
relies on observations of intermediate state, low- disk-fraction black hole binaries. These 
systems exhibit un usual temperature-luminosity relations akin to the ones seen in LMC X-l 



( JDunn et al. 



20111 ). One can confirm whether the same mechanism is at work in both 
LMC X-l and an intermediate state system by measuring the variability (or lack thereof) 
of the latter's total photon flux. Stable photon fluxes would support the proposition that 
these systems also contain stable, but variably-occulted, inner disks. 

LMC X-l has provided specific evidence that stellar mass black hole binaries may have 
x-ray spectral features which are unique to wind accretion. In a more general sense, the fact 
that LMC X-l shares certain x-ray spectral properties with other types of black hole binary 
systems motivates further investigation of its seemingly unorthodox accretion mechanism's 
full scope and utility. 

This work was supported by NASA grant NNX09AC86G. The authors thank the 
anonymous referree for his or her constructive and insightful suggestions. L.R. thanks E. 
Ramirez-Ruiz, R. Strickler, and J. Naiman for productive discussions. 



-25- 



Observation 
ID 


UT (yyyy- 
mm-dd) 


Start Date 
(MJD) 


Total Expo- 
sure (s) 


Sub- 
Observations 


20188-01-02-00 


1996-12-30 


50447.438 


9919 


4 


20188-01-03-00 


1997-01-18 


50466.347 


9941 


3 


20188-01-05-00 


1997-03-09 


50516.354 


10189 


4 


20188-01-06-00 


1997-03-21 


50528.004 


10625 


4 


20188-01-07-00 


1997-04-16 


50554.193 


11484 


3 


20188-01-14-00 


1997-09-09 


50700.695 


10174 


4 


20188-01-18-00 


1997-10-10 


50731.611 


11436 


3 


30087-01-02-00 


1998-01-25 


50838.667 


10017 


3 


30087-01-03-00 


1998-02-20 


50864.370 


9948 


4 


30087-01-04-00 


1998-03-12 


50884.495 


10007 


7 


30087-01-06-00 


1998-05-06 


50939.149 


11396 


4 


30087-01-07-00 


1998-05-28 


50961.101 


5966 


2 


30087-01-08-00 


1998-06-28 


50992.097 


9921 


4 


30087-01-09-00 


1998-07-19 


51013.909 


10094 


4 


80118-01-06-02 


2004-01-11 


53015.340 


10908 


3 


80118-01-07-00 


2004-01-12 


53016.258 


10788 


3 


80118-01-08-00 


2004-01-10 


53014.236 


6608 


2 



Table 2:: Summary of the deep archival RXTE observations used to constrain the power 
law index, and to find upper limits on the apparent timescale of disk luminosity evolution. 
The former used only the parts of the spectra above 12 keV, while the latter relied only on 
regions of the spectra below 12 keV. 



-26- 







LMC XI 


LMC X -3 


Inputs (frozen) 


nH (cm 2 ) 


7.2 x 10 21 (1) 


3.2 x 10 20 (1) 


r 


2.68 


2.34 (2) 


Scattering Flag 


1 


1 


Outputs (average values) 


kT e (keV) 


0.90 ±0.06 


1.07 ±0.18 


Jsc 


0.13 ±0.07 


0.08 ±0.17 


X 2 /d.o.f. 


0.79 ±0.24 


0.79 ±0.77 



Table 3:: Input parameters and results for the phabs x (simpl x diskbb) fits to LMC X- 
1 and LMC X-3 data, with citations in parenthesis. Uncertainties in the outputs correspond 
to single standard deviations of the values about their average. Fixi ng simpl' s scatterin g 



fla g at a value o f 1 for ces the fits to account for only up-scattering. (1) 



Nowak et al, 



(2) 



Smith et al 



(120071) 
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Fig. 5. — : The innermost disk temperature plotted as a function of observation date for 
LMC X-3 (top) and LMC X-1 (bottom). 



-27- 



c 

o 

c 

T3 
Ld 



o 

c 







0.80 0.85 0.90 0.95 1.00 1.05 
Disk Temperature (keV) 



1.10 



Fig. 6. — : LMC X-l inner disk temperature vs. luminosity over the course of two repre- 
sentative archival observations. The top panel shows observation 20188-01-05-00, while the 
lower panel shows observation 30087-01-04-00. Black points, representing subsections of the 
archival observation, are numbered in chronological order. Gray points represent the data 
from our own observing campaign, and provide a basis for comparison. 
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Fig. 7. — : From top to bottom: disk luminosity, coronal luminosity, total photon emission 
rate, and total luminosity of LMC X-lvs. observation date. Luminosities are shown in terms 
of the source's Eddington luminosity. 
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Fig. 8. — : The fraction of total disk seed photons that undergo Compton up-scattering 
versus the inner disk temperature. 
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Fig. 9. — : The obscured inner disk model, fitted to all of the chopped archival observations 
containing more than three sub-intervals. Observation numbers are listed in the upper left 
corner of each plot. Black points represent the chopped observations, while gray points 
represent our data. The solid black line plots the obscured inner disk model, which has been 
normalized to fit the chopped observations. 
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Fig. 10. — : A cartoon of the accretion geometry of wind accretion systems. Incoming mate- 
rial from the companion's wind (dashed gray lines) must dissipate its angular momentum in 
a shock front behind the compact object, a.k.a. the 'post shock', before it can be accreted 
(thick arrow). 
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